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We predict theoretically and demonstrate experimentally an ellipticity-dependent nonlinear 
magneto-optic rotation of elliptically-polarized light propagating in a medium with atomic coher- 
ence. We show that this effect results from hexadecapole and higher order moments of the atomic 
coherence, and is associated with an enhancement of Kerr and higher orders nonlinearities accom- 
panied by suppression of the other linear and nonlinear susceptibility terms of the medium. These 
nonlinearities might be useful for quantum signal processing. In particular, we report an observa- 
tion of an enhancement of the polarization rotation of elliptically polarized light resonant with the 
5S 1/2 F = 2 ^ 5P 1/2 F = 1 transition of 87 Rb. 
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I. INTRODUCTION 

Nonlinear magneto-optic rotation (NMOR) of the po- 
larization plane of light resonant with atomic transitions 
is attracting increasing attention jl], || . Ultranarrow (up 
to 1 Hz H, ||) spectral features accompanied by strong 
polarization rotation observed in NMOR experiments are 
used (or proposed to be used) in sensitive magnetome- 
try |6[ H, in time-reversal-invariance violation experi- 
ments |§, K j m measurements of the electron dipole 
moment [0, 12], and in measurements of various atomic 
constants |1S | . Extremely slow propagation of light has 
also been observed in NMOR in hot rubidium vapor jl4| . 

The most accurate description of the properties of 
NMOR signals is obtained from an analysis of density 
matrix equations for the atomic polarizations and popu- 
lations along with Maxwell equations describing propaga- 
tion of the electromagnetic fields in the atomic medium. 
The exact solution of this problem, however, is very com- 
plicated, and for most cases may be obtained only nu- 
merically. The problem should be somehow simplified to 
obtain analytical results. 

The traditional approach to solution of the problem 
is based on the approximation of weak electromagnetic 
fields and low atomic vapor densities |l5|, [l?], [l8] , con- 
ditions found in early experiments involving incoherent 
radiation from atomic discharge lamps. In this case one 
can use perturbation theory, and the atomic susceptibil- 
ity may be decomposed in a series of the electromagnetic 
fields involved. Magnetic field dependent terms of the 
susceptibility decomposition which are nonlinear in the 
electromagnetic fields are responsible for NMOR. It can 
be demonstrated that only two-photon processes are im- 
portant in this approximation, and therefore complicated 
multilevel systems may be reduced to systems with small 
level number (such as A, V, or A-schemes) lll| pp|. 



'Electronic address: i.novikova@osa.org 



In this approximation, NMOR is a consequence of low 
frequency ground-state coherence formed by two-photon 
processes between Zeeman sublevels with difference in 
magnetic quantum numbers equal to Am = ±2. 

In some cases it is convenient to describe the atom- 
light interaction from the point of view of light-induced 
multipolc moments of the atomic electron distribution. 
Conventionally this is done in terms of an irreducible 
tensor representation of the density matrix |2l], |22|, p3[ . 
In this case, the ground-state coherence is equivalent to 
the quadrupole moment, or alignment. It has been sug- 
gested that NMOR is a consequence of the alignment 
to orientation conversion p4j , where the orientation is 
equivalent to the population difference between nearest 
Zeeman sublevels with Am = ±2. 

The simplified theoretical approaches used for weak 
electromagnetic fields generally fail for strong ones. The 
question that arises here is whether or not the inter- 
action with strong fields bring new physics, e.g. if the 
higher order atomic coherences influence NMOR. Alkali 
atoms have a level structure which allows for a forma- 
tion of the coherent superposition of the magnetic sub- 
levels with Am = ±4 (hexadecapole moment in the mul- 
tipole decomposition of the interaction process) and even 
higher. Such coherences should be excited by multipho- 
ton processes that include four or more photons. Gawlik 
et al. f25f| observed strong narrow features in a forward 
scattering experiment with free sodium atoms, which 
were attributed to a hexadecapole moment. However, 
subsequent work of Giraud-Cotton et al. jl5j and other 
groups [l7| [H], demonstrated that these features may 
be explained using third-order perturbation theory which 
includes only quadrupole moments. 

There have been a number of publications where ob- 
servation of hexadecapole and higher order moments is 
reported for the case where the magnetic field is perpen- 
dicular to the light propagation direction ]^6[ |27| ]. At 
the same time, the question of their influence on forward 
scattering and NMOR signals in Faraday configuration is 
still open |28| . Generally, the interpretation of the exper- 
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imental results in the case of strong laser fields and large 
multipole moments is very complicated. The high-order 
coherence causes only slight modifications of the rota- 
tion caused by the quadrupole moment, which hinders a 
convincing demonstration of these high-order effects. 

We here solve both analytically and numerically the 
problem of the propagation of strong elliptically polarized 
electro-magnetic fields through resonant atomic media in 
the presence of a magnetic field. We particularly investi- 
gate the properties of the light which interacts with the 
magnetic sublevels in an M-like level configuration and, 
therefore, forms coherences with Am — 4. We demon- 
strate that these coherences is responsible for a new type 
of polarization rotation which depends on both the light 
cllipticity and the applied magnetic field. We observe 
this effect in hot vapor of rubidium atoms. Since such 
rotation does not appear for an isolated A scheme, our 
experiment may be treated as a clear demonstration of 
the hexadecapole moment of atoms. 

Another interesting and important feature of the sys- 
tem under consideration is connected with a large Kerr 
nonlinearity that is associated with NMOR. We analyze 
Kerr nonlinearity in the M level configuration and show 
that the ratio between the nonlinearity and the absorp- 
tion may be large. Moreover, we show that by increasing 
the number of Zeeman sublevels (e.g. by using another Rb 
isotope or different alkali atom with higher ground-state 
angular momentum) it is possible to realize higher orders 
of nonlinearities. Our method of creation of the highly 
nonlinear medium with small absorption has prospects 
in fundamental as well as applied physics. It can be used 
for construction of nonclassical states of light as well as 
coherent processing of quantum information (29]. 

To bridge between this and previous studies we should 
note that NMOR may be attributed to coherent popu- 
lation trapping (CPT) |50| pl"| and electromagnetically 
induced transparency (EIT) |p|. Both EIT and CPT 
are able to suppress linear absorption of resonant mul- 
tilevel media while preserving a high level of nonlinear 
susceptibility |||, p5| ]. Previous theoretical studies 
of coherent media with large optical Kerr nonlinearities 
have described nonlinearities resulting from the effective 
self-action of an electromagnetic field at a single photon 
energy level, such as a photon blockade |3(| [37| |3^, |39[ , 
or an effective interaction between two electromagnetic 
fields due to refractive J|4|, ||^, fl0|, ^IJ and absorp- 
tive Kerr nonlinearities. The absorptive \^ non- 
linearities were studied experimentally for quasiclassical 
cases 43 , 44] . It was shown quite recently, that a similar 
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approach may lead to achievement of even higher orders 
of nonlinearity |45[ . 

A method of producing Kerr nonlinearity with van- 
ishing absorption is based on the coherent properties of 
a three-level A configuration (see Fig. |l]a). In such a 
scheme the effect of EIT can be observed. Two optical 
fields, oti and Qi, resonant with the transitions of the A 
system, propagates through the medium without absorp- 
tion. However, because an ideal EIT medium does not 



FIG. 1: Energy level schemes for (a) A-system; (b) TV-system; 
(c) M-system. 



interact with the light, it also can not lead to any nonlin- 
ear effects at the point of exact transparency . To get 
a nonlinear interaction in the coherent medium one needs 
to "disturb" the EIT regime by introducing, for example, 
additional off-resonant level(s) (level 0,2 in Fig. |l]b). In 
the following we refer to the resultant level configuration 
an A^-type scheme. Such a scheme has been used in pre- 
vious works {H H, H> HI EH EH • If the disturbance 
of EIT is small, i.e., the detuning A is large, the absorp- 
tion does not increase significantly. At the same time, 
the nonlinearity can be as strong as the nonlinearity in a 
near-resonant two level system. 

This paper is based on the existence of CPT in multi- 
level media. Unlike the early ideas of Kerr nonlinearity 
enhancement, we propose to use not a single A scheme, 
but several coupled A schemes. In particular, we consider 
the M-type configuration as shown in Fig. |l]c. Coherent 
population trapping exists in such a scheme, like in a 
A-type level system. 

By introducing a small detuning, <5, we may disturb 
this CPT and produce a strong nonlinear coupling among 
the electromagnetic fields interacting with the atomic 
system, while having small absorption of the fields ]^6| . 
The dispersion of the M level media and associated group 
velocity of light propagating in the media are intensity 
dependent due to the nonlinearity, as was theoretically 
predicted by A. Greentree et al. ]^J. Finally, in the 
case discussed below, energy levels of the M configura- 
tion correspond to Zeeman sublevels of alkali atoms. The 
multi-photon detuning is introduced by a magnetic field, 
resulting in the intensity dependent polarization rotation. 

We show a simple way to reduce a five-level M config- 
uration to a four-level N configuration, and prove that 
these completely different schemes demonstrate refrac- 
tive nonlinearities of the same magnitude. This is a very 
interesting result, because the nonlinearity of the M con- 
figuration is a consequence of the hexadecapole part of 
atomic coherence, while the nonlinearity in the N con- 
figuration results from quadrupole atomic coherence. 

Our paper is organized as follows. In Sec. II we ana- 
lyze the F = 1 — > F' = atomic transition, demonstrate 
that this transition may be described by a A level config- 
uration, and show that the polarization rotation in the 
case of a A configuration does not depend on the light 
ellipticity. In Sec. Ill we study F — 2 — > F' — 1 atomic 
transitions, show that it consists of A and M schemes, 
investigate properties of the M interaction scheme, and 
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FIG. 2: (a) Schematic interaction of an electromagnetic wave 
with an atomic transition [&), F = 1 — > \a), F' = 0. The 
electromagnetic field is decomposed into two circularly polar- 
ized components having Rabi frequencies S7+ and f2_. (b) 
Simplification of the scheme (a) for the case when there is 
a magnetic field applied parallel to the wave vector of the 
electromagnetic wave. 



show that ellipticity dependent NMOR is possible. Us- 
ing analytical calculations we show that the hexadecapole 
moment plays an important role here. In Sec. IV we ex- 
pand our theory to the case of generalized M energy level 
systems and discuss possibilities of observations of x' 5 ^ 
and higher order nonlinearities. In Sec. V we discuss ap- 
plications of the nonlinearities for quantum information 
processing. The case of Doppler broadened A, M and N 
systems is considered in Sec. VI for the particular case 
of a weak probe field. We present experimental measure- 
ments of the polarization dependent NMOR in hot Rb 87 
and Rb 85 atomic vapors in Sec. VII. Finally, in Sec. VIII, 
we present our conclusions based on these results. 



emerges as a consequence of this interaction. 

In the following, we briefly review the basic proper- 
ties of CPT in A systems and calculate the optical losses 
and the polarization rotation by solving the optical Bloch 
equations for the density matrix elements. Finally, we 
note how the F = 1 — * F' = level configuration can 
be reduced to a A system via proper renormalization of 
decay rates and density matrix. 



A. Coherent population trapping in a A system 

The Hamiltonian for the A system shown in Fig. 
can be written as 



H A = hA\a)(a\-h5\b + )(b+\ + hS\b-)(b- 
+ h(n-\a){b + \+n + \a){b-\+H.c.) 



(1) 



where E+ and E- are the electric field amplitudes of 
two opposite circularly polarized electromagnetic waves, 
fi_ = E-Pab+ffr, £1+ = E + p a b-/Ti are the correspond- 
ing complex Rabi frequencies, p a b+ and p a b- are the 
atomic dipole moments, A is the one-photon detuning of 
the laser frequency from the exact atomic transition, and 
<5 is the shift of the ground-state sublevels resulting, for 
example, from interaction with a magnetic field. 

The eigenvalues of this Hamiltonian Xi (where H\X) = 
h\\X}) may be found from 
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or 



II. ANALYSIS OF NMOR FOR THE CASE OF 
ANF = 1^F'=0 TRANSITION 

A three-level A configuration is the simplest system 
that results in NMOR. This system appears naturally in 
the configuration of Zeeman sublevels of an F = 1 — > 
F' = atomic transition, where F and F' are the to- 
tal angular momenta of the ground and excited atomic 
states, respectively. This scheme can be easily seen if 
the angular momentum quantization axis is chosen along 
the light propagation direction. The effective interaction 
scheme for this case is shown in Fig. ||a. The A configu- 
ration consists of two circularly polarized components of 
the laser field, which create the low-frequency coherence 
between magnetic sublevels m = ±1. Because of the se- 
lection rules, the electromagnetic waves do not interact 
with the sublevel having m = 0. 

For zero magnetic field such a configuration demon- 
strates coherent population trapping. A nonzero mag- 
netic field collinear with the wave vector of the light leads 
to Zeeman shift of magnetic sublevels m = ±1, which 
disturbs CPT and results in an interaction between the 
light and the atoms. The nonlinear polarization rotation 



A 3 + A 2 A + A((5 2 

5(5A + \fl-\ 2 



fi+l 2 ) = o . 



(3) 



In the degenerate case (6 = 0) the eigenvalues and 
corresponding eigenstates are 



X D 
\D) 





n+\b + ) - 



(4) 



Xbx 



\Bl,2) 



= A± 



A 2 
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n* 



\b+) 



(5) 



The state denoted as \D) is called the "dark state" be- 
cause an atom in this state does not interact with the 
light fields and, therefore, does not fluoresce. Atoms in 
the other two states, called "bright states", readily ab- 
sorb light. Therefore, atoms initially prepared in a bright 
state are optically pumped into the dark state after some 
finite time comparable with the lifetime of the excited 
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level | a). Thus, in steady-state, the atomic ensemble does 
not interact with the electro-magnetic fields, which is the 
essence of CPT. The dispersive properties of the atomic 
system in the dark state are governed by the coherence 
between the ground states of the A system. The corre- 
sponding density matrix element may be found from (see 
Ref. @): 



Pb+b- 



(6) 



The true dark state exists only for 6 = 0. As soon as 
the exact resonant conditions are disturbed, the system 
starts interacting with light. However, for small detun- 
ings (v/|ft+| 2 + \fl^\ 2 > \5\, y/\AS\) the disturbance of 
the dark state is small, and most of the atomic popula- 
tion is concentrated in the modified dark state \D). In 
this case the eigenvalue A_d corresponding to this state 
can be found by solving Eq. (||) and keeping only the 
terms linear in 5: 



An = 5 



l^+l 2 



\D) 



,f2+| 2 + |f2_| 2 
Af\\D) + 25 



(7) 
(8) 



where J\T c± 1 + 0(S 2 ) is a normalization constant. From 
Eq. (||) it is obvious that the population of the excited 
level | a) is proportional to S 2 . 



where N is the density of the atoms in the cell, and v 
is a carrier frequency of the electromagnetic wave. Us- 
ing Eqs. (11) and (0) (with H ~ fiXn) we arrive at the 
following propagation equations for the Rabi frequencies 



dz 



\n T \ 2 



la. 



\2\2 



where k is a coupling constant given by 

3 



8tt 



(12) 



(13) 



and A is the wavelength of the light in vacuum. It is 
also useful to rewrite the equation of motion for the field 
amplitudes E±: 



dE v 
— = T 4inhSN-E ±T 



\E T \ 2 



(14) 



Equation ( |l2"| ) is suitable for describing the phase evo- 
lution of the electromagnetic fields. However the decay 
processes responsible for the optical losses cannot be cor- 
rectly included in this method and we need a density 
matrix approach. In the following section we explicitly 
calculate the density matrix elements for the A system to 
verify Eq. (|l2|) and discuss the attenuation of the light. 

C. Density matrix approach 



B. Equations of motion 

It is possible to obtain the equation of motion for 
the electro- mag netic fields, using the method reported 
in Ref. H^j. If we assume a small disturbance of 
CPT, almost all atomic population remains in a dark 
state during the interaction process, and we can rewrite 
the interaction Hamiltonian as 



H~h\ D \D)(D\ 



(9) 



Since now \D)(D\ ~ 1, the atomic degrees of freedom 
may be excluded from the interaction picture and we 
can write H ~ hXjj . The interaction Hamiltonian may 
be rewritten in the Heisenberg picture, so that Q oc a, 
where a is the annihilation operator for the electromag- 
netic field pq |. The quantum mechanical equation for 
the electromagnetic creation and annihilation operators 
may be presented in the following form: 



da 



i dH 



(10) 



The propagation equation for the electromagnetic field 
amplitude E can be obtained from Eq. (^) as a quasi- 
classical analogue of Eq. ( |l(i| ) Q : 



dE „ v dH 
— =2mN- — 

oz c oE* 



(11) 



In order to discuss a realistic model of the atom-field 
interaction in an atomic cell we need to include atomic 
level decay rates (Fig. ||b). We introduce the decay rate 
7o outside of the system that is inversely proportional to 
the finite interaction time of the atoms and electromag- 
netic field. This decay represent the atoms leaving the 
interaction region. Another term that describes decay 
to outside levels, j r , describes population pumping into 
states that do not interact with the fields, for example, 
the state with zero magnetic moment (m — in Fig. m) . 
The natural decay rate from level \a) to levels |6+) or 
|6_) is denoted as 7,-. 

We also need to take into account the atoms entering 
the laser beam. To do that, we include incoherent pump- 
ing to all Zeeman sublevels from outside of the system, 
which means that atoms that enter the interaction region 
have equal populations of the ground state sublevels and 
no coherence between them. The value of the incoherent 
pumping rate is chosen to be 70/2 to keep the sum of 
level populations equal to unity in the case of j r = 0. 
When 7,, ^ 0, the sum of the populations is less then 
unity because of the optical pumping, i.e., 



Pan + Pb+b+ + Pb-b- 



1 > 
= 1 Paa 

70 



(15) 



The time-evolution equations for the density matrix 
elements pij for the A system can be obtained from the 
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Liouville equation: 



P=-^[H^p]- l -{T lP } + R (16) 

where p = H A is given by Eq. (@), T is the 

matrix describing the decays in the system, and R is 
the matrix of incoherent pumping to the ground state 
sublevels. Then the equations for the atomic populations 
are: 



7o 

Pb- b - = —- 70Ph- b- + IrPa a 

+i(Q* + p ab - - c.c.) 
7o 

Pb+ 5+ = y~ 7oPfc+ b+ + IrPa a 

+i(Q*_p ab+ - c.c.) . 
Analogously, for the polarizations we have 

Pab± = —F ab ±p ab ±+itl T (pb±b±—Paa) 

+iQ±p bT b ± 

Pb-b+ = —Fb-b+Pb-b+ + iQ*+Pab + 
-iQ-Ph-a 



where 



r a6 ± = j + i(A±S) 
r&- b+ = 7o + 2z<5 



(17) 
(18) 

(19) 
(20) 



(21) 
(22) 



with 7 = 7 r + 70 + 7r/2. 

In the steady state case, we can solve Eqs. ( |l9| ) and 
( ^oj ) in terms of the atomic populations: 



Pb- b+ 



Pa b± 



r>b- 



rib-\ 



r&-5+ + 



in. 



in. 



(23) 



b- a 



rib± a I r^zphi 



|n T ] 2 
rb=p Q 



Tlb^a- 



ab± 



Tb T b± 



r„ 



|n T | 2 



(24) 



where n b ± a = p b ± b ± — p aa . Inserting these expressions 
into Eqs. (0) and ( [l8| ) and using the condition given in 
Eq. (15|) we can derive linear equations for the atomic 
populations. In the general case, however, their solution 
is very cumbersome. 

Let us consider the case of a strong electro-magnetic 
field, such that |f2| 2 /7o7 ^S> 1. We also assume that 
|<5|,7o <C 7, |n|, and A = 0. In the zeroth approxima- 
tion the atomic populations are determined by Eq. (Q): 



(o) 
Pb± b± 

JO) 
I a a 



(25) 
(26) 



where \Vl\ 2 = \VL+\ 2 + |ft_| 2 . 

Now we can solve for the polarizations p a b± , keeping 
only the terms linear in S and 70 



Pa b± 



gi(!W±2z^ ± | 2 



(27) 



It is important to note that this expression for the po- 
larization, obtained for an open A system, coincides with 
the analogous expression calculated by Fleischhauer et 
al. H for a closed system, if the ground-state coher- 
ence decay rate and the population exchange rate be- 
tween ground states are the same and equal to 70. This 
proves the equivalence of the open and closed models for 
the description of A schemes, which has been previously 
demonstrated by Lee et al. E{| for the particular case of 
a weak probe field. 

The stationary propagation of two circularly polarized 
components of the laser field through the atomic medium 
is described by Maxwell-Bloch equations for the slowly- 
varying amplitudes and phases: 



dn± 

dz 



in 



(28) 



the 



Note that Eq. (12) can be obtained from Eq. 
limit 70 = 0. 

Separating the real and imaginary parts of Eq. 
and using Q± = \Sl±\ e l ^ ± , one can find the propagation 
equations of the electromagnetic field intensity |f2| 2 and 
the rotation angle of the polarization ellipse <f> — (4'+ ~ 
0_)/2: 



d\n\ : 

dz 



dz 



-K70 
2kS 



(29) 
(30) 



After integration, the following expressions for the light 
transmission J out and the polarization rotation angle cf> 
are obtained: 



-^out — lit 



1 



|n(o)| 2 



— l n 

70 lout 



(31) 
(32) 



where L is the interaction length. It is important to note 
that the final expressions in Eqs. (31) and j3^ ) include 
only the total laser intensity, not the intensities of the 
individual circular components. This means that both 
transmission and polarization rotation are independent 
of the initial polarization of light |50) . 



D. Normalization conditions for an F — 1 — » F' = 
transition 

The correspondence between the F = 1 — > F' = 
scheme (Figs. |^a) and A scheme (Fig. |^>) can be obtained 
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if we exchange 7 r by 7aa/3, where j aa is the decay rate 
of the excited state to the ground state. The decay rate 
7,. should be presented as "f r = 7 aa /3 + 7 aa , where 7 aa 
stands for the decay of the excited state outside of the 
system in Fig. ^|a. 

We assume that the incoherent pumping rate into each 
Zeeman ground state is equal to 7o/3, to keep the nor- 
malization condition similar to Eq. (|15|): 

7a a 

Paa + P+1, + 1 + P-l.-l + P0,0 = 1 —Paa (33) 

7o 

where p aa is the population of the excited state and pa is 
the population of the ith magnetic sublevel of the ground 
state in the system depicted in Fig. ||a. 

Keeping in mind that the population of the m = 
state is determined by the decay rate of excited state \a) 
and by the decay outside of the system we write the rate 
equation 




FIG. 3: a) Energy level scheme for 87 Rb atoms. This scheme 
may be decomposed into a superposition of b) A-system and 
c) M-system. Transition probabilities are shown for each in- 
dividual transition. 



70 , laa - 

PQ,0 = "g 70P0 : + —Paa 



and solve it in the steady state 



Po,o 



1 | laa - 
o T „ Paa 

3 37o 



(34) 



(35) 



Let us assume that p aa = £p aa , p+x t+ i = £,Pb+b+, and 
= £,pb-b-- The normalization parameter £ can be 
found by substituting Eq. (|35| ) into Eq. j3^) , and com- 
paring the normalization conditions Eqs. (|15|) and (pq): 



(36) 



Therefore, we can derive density matrix elements for the 
F = 1 — » F' = level scheme shown in Fig. ||a by simple 
multiplication of the elements of the density matrix for 
the A scheme by the scaling factor £. 



M scheme and a A scheme is that the higher order coher- 
ence (Am = 4) becomes important. Since the A system 
had been studied in the previous section, we primarily 
concentrate on M scheme here. 

The M scheme is described by a set of twelve den- 
sity matrix equations. The only straightforward way to 
solve this system is with numerical methods. However, 
if we study the atomic interactions with weak magnetic 
fields, the decay processes and polarization rotation pro- 
cesses are independent, as we saw for the A configuration. 
Thus, the polarization rotation may be found in analyti- 
cal form under the condition of zero relaxation using the 
Hamiltonian diagonalization procedure as presented for 
the A system. The modified Schrodinger equation model 
is suited for this as well. The optical losses may be found 
separately by considering the optical pumping into the 
dark state with zero magnetic field. 



III. ANALYSIS OF NMOR FOR THE CASE OF 
AN F = 2 -> F 1 = 1 TRANSITION 

For atomic ground atomic state angular momentum 
higher than F = 1 it is possible to create more than one 
A link between magnetic sublevels. This is equivalent to 
the creation of coherent atomic states characterized by 
higher angular momenta, which may drastically change 
the interaction of such a medium with the electromag- 
netic field. 

Let us concentrate first on F — 2 — > F 1 = 1 transitions, 
which occur in the 87 Rb Di line. The case of higher angu- 
lar momenta is discussed in the next section. Interaction 
of elliptically polarized light with the F = 2 — * F' = 1 
transition may be decomposed into a A scheme with 
m = — 1 <-> m' = <-> m = +1, and an M scheme 
m = —2 <->m' = — l<->m = 0<->m' = +1 <-> m = +2, 
as shown in Fig. Km. The main difference of between an 



Coherent population trapping in an M level 
scheme 



It has been shown that the dark state exists even for 
atoms with complicated Zeeman substructure interacting 
with elliptically polarized light [|lj 



55, 56 



_ 3 B Ei . 

Here we recall the analytical expressions for this dark 
state and the corresponding eigenvalues. Using an ef- 
fective interaction Hamiltonian, we derive propagation 
equations for the electromagnetic fields. We restrict our 
consideration to the case relevant to the M configura- 
tion consisting of Zeeman energy sublevels in the mag- 
netic field. That is, we assume that the atomic tran- 
sition frequencies are such that w a -bo = f-^a+bo = 
Lu a -b- — w — 25, and u> a +b+ = lo + 25, where the detun- 
ing S is due to a Zeeman shift, and the laser frequency v 
is resonant with the atomic transition. The interaction 



7 



Hamiltonians for M systems is 



H M = -2hd\b+)(b + \ 

+ h(ni-\ a+ )(b^ 

+ n 1+ \a+){b \-\ 



+ 2HS\b-)(b-\ 
\ + n 2 -\a-)(b \ 
n 2+ \a-)(b-\ + 



H.c. 



(37) 



where 0i_ = E-p a+b+ /h, = E + p a + b0 /h, fi 2 _ = 
E_p a _ b0 /h, n 2 + = E + p a _ b _/h (see Fig. |c). 

As in the A system, the eigenvalues of the interaction 
Hamiltonian can be determined from: 



25 



-A 


o* 








h+ 


-A 










n* 2 _ 


-A 


o* 

"1+ 








fii+ 


-A 


















Ql_ 

-A — 25 



(38) 



Also, as in the A scheme, the eigenvalue A = and cor- 
responding "dark state" exists only for 5 = 0: 



\D) 



n 1+ n 2 +\b+) - fti_f7 2+ |6 ) + fii_o 2 _|&_) 



2+1 



|o 1+ | 2 |ft 2+ | 2 + |f>^| 2 |fv 



( 39 ) 

It is worth noting that, similar to the A system, the non- 
vanishing low frequency coherences p b +w and p b Q b - are 
important here. The major difference in the dispersive 
properties of the M and A schemes arises from the exis- 
tence of the four-photon coherence p b + b - : 



Pb+b- 



r2 1 _r^2_^i+^ 2 + 



i+I 



2+1 



For small 5 we can again find the eigenvalue for the quasi- 
dark state, taking into account only the linear terms in 
5: 



A 



M 



25 t 



|fi 2 -| 2 |tti-| 



l^+H^i+l 2 



,f! 2+ | 2 |r! 1+ | 2 + |fi 2 -| 2 |fti-| 2 + |n 2+ | 2 |n!_| 2 • 

. (41) 

Using Eq. (XViY) we derive equations of motion for the 
fields. For example, 



dfl 2 



- = 2inSQ 2 - 

2|n 1+ i 2 ir! 1 _ 



3 a-b0 



J 2+| 



(42) 



I0 2 



[Hi-I 



|f2i-| 2 |« 2 +| 



|«l-| 2 |^2-| 2 ) S 



where k is the coupling constant with respect to the 
transition as a whole (i.e., j r in Eq. ( |l3| ) is now the 
total natural decay rate of the excited state), and p = 
(4z/ 3 7 r /(3?ic 3 )) 1 / 2 is the dipole moment of the transition. 

The calculations can be considerably simplified if the 
numerical values of the transition probabilities are used. 
Let us now consider the particular case of the M part 
of the F = 2 — > F = 1 transition. According to 
the transition probabilities, shown in Fig. pic, we get 
N2+I7N2-I 2 - 6|£ + | 2 /|£L| 2 and |0 1+ |7|fi 1 _| 2 = 
|-E + | 2 /6|_E_| 2 . The interaction Hamiltonian (Hm — 



TiXm) for the elliptically polarized laser field can therefore 
be rewritten as 



H M ~ 2h5- 



|£-| 4 -|£+| 4 



\E + \± + \E^ + 6\E + \2\E_ 



(43) 



and therefore 

dE± 
dz 



T8iirh5N-E± x 



, 3(\E + \ 4 + \E^) + 2\E + \ 2 \E^ 

(\E + \* + \E^ + 6\E + \i\E^y 



.(44) 



In what follows we derive the same equation using the 
more rigorous modified Schrodinger formalism [B7|. 



B. Solution based on the modified Schrodinger 
equations 

The interaction described above of the four electromag- 
netic fields with the M energy level configuration may be 
also studied using Schrodinger equations. This approach 
enables us to find exact expressions for all the atomic ob- 
servables when we can ignore spontaneous emission. The 
state vector of the atom can be written as: 



I*) 



a+e 
bo\b ) 



l \a + )+a-e- ivt \a-) + 



(45) 



b+\b + )+b_\b_). 
Solving the Schrodinger equation 

I*) = ™#|*> 

for the interaction Hamiltonian Eq. fl37|), we obtain the 
following equations of motion for the slowly- varying state 
amplitudes: 

(46) 
(47) 
(48) 
(49) 
(50) 

In the steady state regime, this system has a nontrivial 
solution only for 5 = 0. The solutions for nonzero detun- 
ings correspond to zero amplitudes for all parameters. 
Thus, to sustain steady state in the open system, exter- 
nal pumping is necessary. For a small splitting between 
ground state levels H5 <C kT, where T is the temperature 
of the vapor, we assume that in thermal equilibrium, i.e., 
in the absence of all fields, all lower states \b±) and 1 60) 
are equally populated. And, therefore, within the open- 
system approach, we assume that the atoms are pumped 
into states |&+), \b-), or \bo) with equal probability from 
outside of the system. The corresponding rate can be 
determined by the requirement that the total probability 
to find an atom in any of the states is unity. 



«+ 


= i^i+bo 


+ jf2x_ 






= iQ 2+ b- 


+ z£l 2 _ 


-b 




= 2i5b+ + 


iQi_a 


+ 


b- 


= -2i5b_ 


+ iVL* 2 ^ 




k 


= iQ* + a+ 


+ iVl* 2 - 
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Unlike the density matrix approach, a straightfor- 
ward introduction of incoherent pumping into the ground 
states of the system is impossible. It was shown by Fleis- 
chhauer j5?j] in an elegant way that the effective density 
matrix equations for open systems with injection rates 
into states and decays out of states can be written in 
terms of stochastic complex state amplitudes. 

Let us consider an effective density matrix equation 
for an atomic ensemble undergoing a unitary interaction 
with some external fields or potentials. In addition, de- 
cay out of atomic states |j) is taken into account with 
rates jj. Also injection into certain states is considered 
with injection rates Rij. In our case the injection occurs 
only into energy eigenstates of the atoms or incoherent 
mixtures of them, so only diagonal elements of the matrix 
Rij are nonzero. If injection in the coherent superposi- 
tion states is considered, non-diagonal elements are also 
required to be taken into account. 

An effective density matrix equation has the following 
structure: 



Pij(t) — Ri 



7i + Ij 



P33 



(51) 



where 7^ are decay rates out of the system, which can in 
general be different for individual states. Generally, the 
pump rates Rij are time dependent, but for the sake of 
simplicity we assume in the following that the rates Rij 
are constant. 

Density matrix elements may be represented in terms 
of state amplitudes pji = c*Cj . In order to put the pump 
term Rij in a similar form, we introduce a formal Gaus- 
sian stochastic variable with the following properties: 



(n) = 
(r irj ) = 

(r* rj ) = R,j 



(52) 
(53) 
(54) 



This yields a set of amplitude equations with stochastic 
pump terms: 



2 l + n 13 3 



(55) 



Since the amplitude equations are linear, their solution 
will be a linear functional of the stochastic pump rates 
j-j. Thus the averaging of bilinear quantities such as c*Cj 
required to obtain the density matrix elements can eas- 
ily be performed. Generally, solution Cj of Eq. ( p5| ) no 
longer makes sense as the amplitude for the atomic wave 
function. It only determines density matrix elements of 
the system. 



To apply the above technique to our problem, we 
rewrite Eqs. (|46|)-([5"o|) (with time derivatives set equal 
to zero) 



iQi+bo + zOi_5 + = 
zf2 2+ 6_ + iQ,2~b$ = 
2i5b + + if2*_a_|- = ir + 



-2i5b- + iO,2+ a — = * r - 



(56) 
(57) 
(58) 
(59) 

jfij + a + + iQ^a- = irQ (60) 
where the stochastic "pumping" is introduced 

<r±) = (r ) = 

(r±r T ) = (r±r Q ) = 

(r±rzf) = (r*±r ) = 

(r±r ± ) = (r*r ) = r 2 . 



Solving Eqs. (|56|)-(|60D with respect to 01, a 2 , b±, and 
60 we get 

-oq— — , o- = -b - 



bn = 



r + |n 2+ | 2 fii_oj, + r_|ni_| 2 n2+n5_ - r |Oi-| 2 |n 2 +| 2 



25 (\n 1+ \ 2 \n 2+ 



r \n 



a+ = — 



|n 1+ |2|n 2+ |2-|n 1 _p|n 2 _|2 

r + \n 2 -\ 2 fli- + r_fi$_f2i + f2 2 + - r |f2 2+ | 2 f2i- 



(61) 



Utilizing the normalization condition 

(a*_a-) + (a* + a+) + (b* + b+) + (btb-) + (b* Q b ) = 1 (62) 
we get 

r = 2d (\n 1+ \ 2 \n 2+ \ 2 - |Q!_| 2 |0 2 _| 2 ) / (63) 

{(|n x _| 2 |n 2+ | 2 + |f> 1+ | 2 |f> 2+ | 2 + \n^\ 2 \n 2 ^\ 2 ) 2 + 

r 'i+i 
2(\n 1+ \ 4 \n 2 



45 2 [\n 1+ \ 2 \n^\ 2 (\n 1+ \ 



'2+| 



1/2 



Nl-I 4 |tt2-| 2 )]} 



Using Eq. ( J63| ) we arrive at the complete solution of the 
problem which takes into account all orders in 6. For 
<5 = the system is in a dark state and the density ma- 
trix elements correspond to the elements generated by 
Eq. (39). For a nonzero small two-photon detuning the 
populations and coherences for the ground state stay ap- 
proximately unchanged. The solution for the populations 
of the excited states are 
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Pa+a+ 



4(5 



3) n 1+ | a |ni_| a |n 2 



|fil + | 4 |«2' 12 



'2+| 



|Ol-| 4 |fl 2 -| 2 



|^1-| 2 |^ 
|2|o. 121 



|0 1+ | 2 |^ 2H 



= 4(52 |fi 2+ | 2 |0 2 _hfi 1+ |^ + |f7 2+ | 4 |0 1+ 



|tti-| 2 |f! 2 -| 2 ) 2 

2 + |r» 2 _| 4 |fi 1 _| 



(|n!_|2|n 2+ |2 + |fi 1+ |2|fi 2+ |2 + |r !l _|2|r! 2 _|2)2 
I 



(64) 
(65) 



and for the atomic polarizations are 



Pa-bO 



Pa+bO 



Pa-b- = 



Pa+b+ — 



25n 2 -{2\n 1+ 


l 2 Ni- 


| 2 |^2+ 




-|Oa+ 


l a |Ox- 


I 4 ) 


(l«i-l 2 l«2+| 2 - 


f |fii+ 


I 2 |^4 


■I 2 h 


-|ni_ 


I 2 |«2- 
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25n 1+ (2\n 1 - 


I 2 |^2+ 


I 2 |^2- 


I 2 ^ 




| 2 |^2+ 


I 4 ) 


(I«l-I 2 |«2 + | 2 - 


f |fil+ 


PI^H 
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-Ifll- 


I 2 |«2- 


I 2 ) 2 


2sn 2+ (2 \n 1+ 


| 2 |^_ 


I 2 |«2- 


l 2 ^ 


-|n a - 


I a |«l- 


I 4 ) 



(I^PI^+P + |fi 1+ |2|fi 2+ |2 + |fi 1 _|2|r! 2 _|2)2 

2^ 1 _(2|^ 1+ | 2 |n 2+ | 2 |» 2 _| 2 + |» 1+ | 2 |» 2+ | 4 ) 

(|fii_| 2 |n 2+ | 2 + |tt 1+ | 2 |tt 2+ | 2 + i^-pi^-i 2 ) 2 ' 



(66) 
(67) 
(68) 
(69) 



r 



Here we kept only the lowest order terms in 5. In the ex- 
pressions for the atomic polarizations, the first term, con- 
taining the amplitude of all four optical fields (for exam- 
ple, f2i + |£!i_| 2 |fi 2+ | 2 |f2 2 _| 2 in the equation for p a +bo), 
is due to the four-photon coherence (hexadecapole mo- 
ment), whereas the second term represents the effect of 
optical pumping. 

The propagation equation for the fields is 



dz 



2-kv 



, r Pij 



(70) 



where the indexes ij show that the values are related to 
the same transition \i) — ► \j). It is easy to see, for ex- 
ample, that the matrix element in Eq. (pq) results in the 
propagation equation in Eq. (f42|). The two approaches 
are therefore equivalent. The equation of motion for the 
circularly polarized electromagnetic fields in E± are given 
by the following expressions: 

dE + ,2irv 

= l ^ {Pa-b-Pa-b- + Pa+bQPa+bO) (71) 

oz c 

dE_ 2tw a . , , . 

-T. = t N (p a +b+pa+b+ + Pa-bOPa-bO) ■ (72) 

oz c 

Substitutin 
Eqs. 



the expressions for atomic polarizations 
and using the proper dipole moments for 
each transition (for the 87 Rb they are equal 1/2 for 
\b±) -> |a±>. and 1/12 for |6 ) -> |o±) (Fig. |c)), we 
obtain Eqs. (Q. 

So far we have made no assumption concerning the 
losses in the system. Generally, this requires solving the 
Bloch equations for the atomic populations and polariza- 
tions as was done for the A system. For the M scheme, 



however, this process is rather involved even for the de- 
generate system (5 — 0). Since the dark state exists 
for any value of Rabi frequency fiy, it is always possi- 
ble to transform the basis of the atomic states so that 
there is one atomic level uncoupled from the laser field. 
The M system can be represented as two independent 
open two-level systems, connected only via relaxation 
processes |5^] . The absorption in this systems has sim- 
ilar properties compared to those of the A system: it is 
proportional to decay rate 70 and inversely proportional 
to the light intensity. The exact analytical expression for 
this absorption is rather lengthy and we do not present 
it here. 



C. Polarization rotation for an F — 2 — > F' = 1 
transition 

To describe the polarization rotation on the F = 2 — > 
F' = 1 transition we write the interaction Hamiltonian 
as a balanced sum of the Hamiltonians for the M and A 
systems, taking into account the branching ratio for the 
atomic transitions 



H 2 ^i =dH A + C, 2 H M = (ih\ A + c 2 h\ 



M 



(73) 



where Ci and C2 are weighting coefficients (Ci + C2 = 1) 
that describe the population redistribution between the 
A and M schemes. Using the numerical simulation of 
this system, we find them to be equal with very good 
accuracy. Using Eq. ( [TT| ) we now derive the equation of 
motion for this system: 
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dE ± 

dz 



^UirhSN-Ei 



\E T \ 2 



E. 



\2\2 



1 + 2(1^4 



r 



|2\2 



3(1^4 



|S_| 4 ) + 2|£ + | 2 | J E;_| 



(1^-4 



|£_| 4 + 6|£; + | 2 |£;_| 2 ) 2 



(74) 



It is interesting to note that for linearly polarized light 
= | f2_ |) the contributions from A and M system 
are identical, and Eq. (Q) coincides with Eq. (0). This 
proves that a single A system may be used for accurate 
description of the dispersive properties of more compli- 
cated level configurations. 

Let us introduce the electromagnetic field ellipticity 
parameter q such that the amplitudes of the circularly po- 
larized components are E± = \E\y/{l ± q) exp(j(/)±)/v / 2- 
Then Eq. (Ff3) transforms to 



dE± 

dz 



T2iirh6N 



v E±(1tq) 
c \E\ 2 



1 



(2 



(75) 



Based of the results of our numerical simulation, we con- 
clude that absorption of light that interacts with the 
F = 2 — * F 1 = 1 transition does not depend on the el- 
lipticity of the light. The light transmission through the 
cell can be described by an equation similar to Eq. (31): 



In 



= Ita 1 - 



2irhj NL v 
|£(0)| 2 c 



(76) 



The rotation angle for the light polarization is then given 
by 



S_ 

To 



1 



(2-<Z 2 



ln- 



(77) 



where I ln and 7 out are the intensities of the electromag- 
netic field at the entrance and exit of the medium. The 
value of polarization rotation increases with the light el- 
lipticity by the factor 



!>M+A 
4>A 



1 



(2- 



(78) 



compared to A system. Therefore NMOR on the F = 
2 — > F' = 1 transition may only be properly described 
by a A configuration for linearly polarized light. The 
difference between the M and A systems results from the 
hexadecapole moment induced in M configuration. 



IV. NMOR IN ATOMS WITH LARGE VALUES 
OF ANGULAR MOMENTUM 

Higher order coherence can be excited among Zeeman 
sublevels of alkali atoms with F > 2. Here we find 
a perturbed dark state for the generalized M scheme 
consisting of an arbitrary number of A-links, using the 
method described above. Then we apply these results 
to evaluate the nonlinear Faraday rotation in the 85 Rb 



F = 3 — > F = 2 transition. We consider the scheme in 
Fig. |^. The interaction Hamiltonian for this scheme is 

71 

H nxA = -nsJ2(n-2k)\b k+1 )(b k+1 \ (79) 

n 

+ (fi fe _|afc)(&fc| + n k +\a k ) (b k +i\ + H.c.) 

fe=l 

Here n is the number of A links, which connects n + 1 
ground-state levels via n excited states. There exists a 
dark state for this system for exact resonance ((5=0): 



n 

B 


-i) fe n>- 


n 

n n '+ i^+i) 


fc=0 


3=1 


l=k+l 




n k 


n 






i 2 n i^+i 2 


\ 


fc=oi=i 


l=k+l 



(80) 



where we use a convention that Yij=i = n?=n+i = 
We deduce the perturbed "dark state" eigenvalue for the 
Hamiltonian using the same procedure as we used before 
in Eq. © 



ixA 



^(2fc-n)[]l^-| 2 II 

i=l j = k+l 



,k=0 



k=0 i=l j=k+l 



(81) 



The equation of motion for the circularly polarized elec- 
tromagnetic fields can be found from Eq. (|Tl"|). As an 
example, let us calculate the interaction Hamiltonian for 
light interacting with the 5S1/2F = 3 — > 5P1/2F' = 2 
transition of 85 Rb (Fig. |J). The circularly polarized com- 
ponents of the resonant electromagnetic field form an M 
scheme and a triple-A scheme. Using the proper values 
of the transition probabilities, shown in the same Figure, 
we derive 



2Ci 



\E\ 4 ~\E 4 



■3E_ 



ioi^pi^-I 2 



+ 



5\E+r 2 \E_ 



5\E + \*\E_\ 2 -\E + f 



15|£-| 2 |£ + | 



15|£ l 4 |£ 4 



\E\ 



(82) 



Here again £1 , 2 are the coefficients reflecting the popula- 
tion distribution between to schemes. By differentiating 
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(a) 



(b) 



\ E_ 



y - 65 



m=l 



FIG. 4: Energy level scheme for Rb atoms. This scheme 
may be decomposed into a superposition of a) M-system and 
b) triple-A system. Transition probabilities are shown for each 
individual transition. 





2n8 



-&1 



FIG. 5: Generalized M interaction scheme. Here Qj_ = 

E-p aib Jh, tt i+ = E+pa^^JH. 



the Hamiltonian it is easy to find the polarization rota- 
tion in the system 



50 
dz 



-6inN- 



hS 

W 



20 



(4-g 2 ) 2 



+ C2 



■ - 6q 2 + 3q 4 
(4-3g 2 ) 2 



It is obvious that both interaction chains contribute to 
the elliptically dependent NMOR. At the same time dif- 
ferent orders of the nonlinear susceptibility are responsi- 
ble for the polarization rotation: if in the case of the M 
scheme it is x^ nonlinearity, for the triple-A scheme it 
is nonlinearity, since there are 7 photons involved in 
the creation of the ground-state coherence. This might 
be the reason why the triple-A scheme shows more en- 
hancement of the polarization rotation for nearly circular 
polarization compared to the rotation of linear polariza- 
tion than does the M scheme (10 vs 20/9 times for the 
F = 3 -> F' = 2 transition). 



APPLICATION FOR QUANTUM 
INFORMATION PROCESSING 



So far we have considered A and M schemes of the 
type described in Figs. ||a and ||c. Here the two-photon 
detunings with respect to states |6+) and |i>_) are equal 
and opposite in sign and all the fields are treated clas- 
sically. This approach is useful for describing NMOR in 
alkali atomic vapors. In general, however, the M system 
may be created by strongly nondegenerate atomic lev- 
els and all four fields connecting corresponding atomic 
transitions may be independent. This case is especially 
interesting if we are going to use the enhanced Kerr non- 
linearity the system provides fl46| . 

In this section we compare the N and M configura- 
tions shown in Figs. [TJd and [l]c. The N system is essen- 
tially a A system with an additional nonrcsonant tran- 
sition. Similarly, the M system in Fig. [j]c is a resonant 
N system with an additional detuned transition. Since 
these systems have potential applications in the field of 
quantum information processing, we discuss them here. 
Some details concerning such systems have been given 
earlier |45|, Q . The systems seems to be completely dif- 
ferent because the all-resonant N configuration demon- 
strates enhanced three-photon absorption, while the all 
resonant M configuration demonstrates complete trans- 
parency. We show here that the performance of these 
schemes as sources of refractive Kerr nonlinearity is very 
similar. 

We assume that the |ai) — ► and \a 2 ) — * \b 2 ) transi- 
tions are induced by quantized fields whereas the transi- 
tions |ai) — * 1 62) and 1 02) — * I&3) are induced by classical 
fields of Rabi frequencies f^i and f^, respectively. 

The Hamiltonians for the N and M schemes in the 
slowly varying amplitude and phase approximations are 

H N = nA\a 2 )(a 2 \ + n(ai\a 1 )(bi\ 

+ Sli\ai){b 2 \+a 2 \a 2 ){b2\+H.c.) (84) 

H M = -M|&3><6 3 | + Maik><&i|+ttiM(&2| 

+ a 2 \a 2 )(b 2 \ +n 2 \a 2 )(b 3 \ + H.c.) (85) 

where H.c. means Hcrmitian conjugate and the relation 
between Rabi frequencies of the probe fields and quantum 
operators describing the corresponding field mode can be 
written as 



'27Tp 2 ^ . 

-di 



hv,. 



(86) 



where pi is the dipole moment of the transition |a,-) — > 
\bi), v.i is the field frequency, V$ is the quantization vol- 
ume of the mode, hi and a\ are the annihilation and 
creation operators. Proceeding along the same lines as 
in Sec. Ill we obtain the effective Hamiltonian for the 
two configurations of the form 



Ti5a\aia\a 2 



(87) 
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where the the coupli ng constant <5 for the two configura- 
tions is of the form |4q| : 



2 9 



and 



3M 



-8- 
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Nil 2 in 2 | 



(89) 



Any system that may be described by the Hamiltonian 
in Eq. ( |§7| ) has a potential application in implementing 
a quantum phase gate. The transformation for a two- 
bit quantum phase gate for the jth and feth qubits is 
given by Q^\aj , fa) = eiq>(ir}6a, :j ,i6p ki i)\aj, fa), where 
| ay) and |/3fc) stand for the basis states |0) or |1) of the 
qubits. Thus the quantum phase gate introduces a phase 
77 only when both the qubits in the input states are 1. A 
representation of the quantum phase gate is given by the 
operator 

Qf = \Oj,o k )(o s ,Ok\ + \oM(OiM 

+ |l„0 fc )(l J; fe |+e"'|l J ,l fc )(l J ,l fc | . (90) 

It is clear that such a phase gate can be realized via 
Hamiltonian H c g with the time-evolution unitary opera- 
tor exp(—iH Q gr /K) and the corresponding phase rj = 5t 
where r is the interaction time. 

The nonlinearities associated with both the present N 
and M schemes correspond to x^ ■ The resonant en- 
hancement of x^ 5 -* and higher order nonlinearities can be 
obtained by adding more A sections to N or M schemes. 
In general, the effective Hamiltonian for x 



(2m-l) 



TT {2m — 1) + 7-rn - t - -t- -+ - 

ii cff = no a[aia2a.2...a' m a ri 



where, for extended N systems, 



Sn = (-1)' 



and, for the extended M system 
(-l) m <5 



if 2 



>M 



Nil 2 |fi 2 
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(91) 



(92) 



(93) 



Such nonlinearities can be used in implementing m-bit 
quantum phase gates that are defined via 



.ou 



exp(ir]S ai ^5 a2t i...S amA )\ai,a2,—a m ) . (94) 

In other words, a phase rj is introduced when all the 
qubits are in state |1). Thus if qubit states |0) and |1) 
arc defined via photon number states, the m-bit quantum 

phase gate is implemented via = exp(— iH m T/h), 

rj = St. Such gates may have important applications in 
quantum computing algorithms such as those related to 
quantum search of unsorted database p9|. 



The important question is how large can the phase shift 
•q be. Our initial estimates indicate that phase shifts as 
large as 3 radians can be obtained for m = 3 via x*- 5 -* non- 
linearities. However, there are problems related to phase 
mismatch between different photons which arise because 
the group velocities can be different for different pulses. 
Such problems can be overcome by methods discussed 



in 1 41 



It is interesting to mention that interaction Hamilto- 
nian for symmetrical M-scheme, given by Eq.(f43|) is iden- 
tical to one for asymmetric M scheme, considered in this 
section - Eq. (j87|) in case when one of the circularly polar- 
ized component is much stronger than the other (nearly 
circularly polarized light). This means that the quantum 
phase gate discussed above, can be potentially created 
even using Zeeman substructure of alkali atoms, resolved 
in magnetic field. Unfortunately, in case of generalized 
M scheme this is not true. 



VI. SUSCEPTIBILITIES FOR 
INHOMOGENEOUSLY BROADENED A, N AND 
M SYSTEMS 

It is important to know what changes are introduced by 
Dopplcr broadening to the systems discussed above. For 
the sake of simplicity we restrict ourselves to asymmetric 
schemes discussed in the previous section. Let us start 
with the Doppler broadened A system shown in Fig. 
This system is widely discussed in the literature 
|6l[ , so we consider only the necessary steps that allow us 
to calculate susceptibility for the Doppler broadened M 
configuration. To sustain EIT in a Doppler-broadened A 
medium the minimum value of the Rabi frequency of the 
coupling field Sli (|Sli| ^> |aii|) should exceed Wd-^/jopy, 
where Wd is the linewidth of the Doppler distribution 

/7o7) f 



(Wd^/70/7 3> ^707) j60fl . Then the population of the 
state 1 61) is almost unity and density matrix element (c.f. 
Eq. |2J) for the probe transition reduces to 




Pabl 



20-1(70 - iS) 



(-/ + i(6 + kv))(j - iS) + inrf 



(95) 



where k is the wave vector of the field, and v is the 
atomic velocity. We simplify the problem by using a 
Lorentzian profile as the velocity distribution function 
f(kv) with full width at half maximum 2W]j such that 
f(kv) = (l/n)W D /[Wp + (kv) 2 ]. Integrating over the 
Dopplcr distribution we get 



(Pa 



bl v 



10-1(70 - iS) 



( 1 + W D -iS)( l0 ~i6) + \n 1 \ 2 
iai 

7 + Wd + i\fli\ 2 /S ■ 



(96) 



This result was evaluated using the contour integration in 
the complex plane which contains one pole in the lower 
half, (kv)i = —iWo- Let us consider the M scheme 
shown in Fig. |T]c ^> \atj\). The susceptibility for 
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the field a 2 may be obtained similarly to the A scheme. 
The population of level \b 2 ) is equal to, approximately, 

lailVinil 8 . 

The nonlinear interaction appears as the result of the 
refraction and absorption of the second probe field a 2 , 
coupled to the second drive field ^2, that create a A 
system. Therefore, we get the susceptibility 



XM 



72(70 ^ iS) 



( l0 - iS)W d + \n 2 \ 2 



(97) 



where N is the atomic density, 72 is the decay rate of the 
level |<i2), A q2 is the vacuum wavelength of the field a 2 . 

Finally, let us consider the N level configuration shown 
in Fig. [fjc. If the condition A ^> 72 is fulfilled, then 
the population of level 1 62) is equal to approximately 
|ai| 2 /|i7i| 2 . The nonlinear interaction appears as the re- 
sult of the refraction and absorption of the second probe 
field «2, far detuned from the corresponding atomic tran- 
sition. For the corresponding two-level system we derive 



Pa2b2 



ia 2 



7 + i(A + feu) l^l 2 " 
The corresponding susceptibility for the field a 2 is 



72 



«1 



(98) 



(99) 



The nonlinear phase shift may be increased, formally, 
by increasing the atomic density or interaction length. 
This is impossible to implement practically because of 
the absorption of the medium. Therefore, to compare 
the nonlinear performance of different nonlinear systems 
one needs to compare the ratio of their refractive nonlin- 
earities and corresponding residual absorption, linear as 
well as nonlinear. The effective ratio between absorption 
and nonlinearity for the Doppler broadened N scheme 
( |99| ) is Wd/A. It is easy to see that ( p^ ) and ( |97| ) are in- 
terchangeble if 70 — *• 0, and A <-> 5/\Ch, | 2 . Therefore, the 
M and N schemes are equivalent in sense of the effective 
Kerr nonlinearity they produce. 



V4 

for ellipticity 
leasurements 



FIG. 6: Schematic of the experimental setup. 



magnetic shield to minimize the influence of the labora- 
tory magnetic field. The atomic density is controlled by a 
heating element placed between the two shielding layers. 
The longitudinal magnetic field is created by a solenoid 
mounted inside the inner magnetic shield. 

To measure the transmitted laser power and the polar- 
ization rotation angle a polarization beam splitter (PBS) 
is placed after the atomic cell. The signals from the two 
PBS channels Si, 2 are collected while the axis of the PBS 
is tilted at 45° degrees with respect to the main axis of 
the initial polarization ellipse. In this configuration the 
transmitted light power is proportional to the sum of the 
two signals Si + S 2 and the polarization rotation angle 4> 
is given by: 



1 . Si — S2 

— arcsm ; 

2 (Si+S 2 )cos2e 



(100) 



It is also possible to detect the ellipticity of the outgoing 
laser beam by placing another quarter waveplate after 
the cell and before the PBS. When the fast waveplate 
axis is aligned with the PBS axis and makes 45° with the 
initial polarization direction, the ellipticity e of the beam 
can be found similarly to the rotation angle: 



1 . S1-S2 
— arcsm -= — 

2 Si + S 2 



B. The experiments with Rb vapor 



(101) 



VII. EXPERIMENTAL STUDY OF NMOR 
WITH ELLIPTIC ALLY POLARIZED LIGHT IN 
RB VAPOR 

A. Experimental setup 

The experimental setup is shown schematically in 
Fig. We use an external cavity diode laser (ECDL) 
tuned in the vicinity of the D\ line of 87 Rb (A = 795 nm). 
The initial linear polarization is produced by a high- 
quality polarizer Pi; the initial ellipticity of the beam 
e |i2| is then controlled by a quarter wave-plate placed 
after the polarizer. Maximum laser power delivered to 
the atomic cell is P m ax — 2 mW. A cylindrical glass cell 
of length 50 mm and diameter 25 mm is filled with iso- 
topically enhanced 87 Rb. It is placed inside a two-layer 



There are two factors contributing to the rotation of 
the elliptical polarization of light propagating through 
the Rb vapor: the nonlinear Faraday rotation, caused 
by the shifts of the magnetic sublevels in an external 
magnetic field, and the self-rotation caused by the ac- 
Stark shifts due to the off-resonant interaction of the 
electro-magnetic field with far-detuned levels [^3[ ^4|, . 
Since the latter effect does not depend on the magnetic 
field, we eliminate it from the experimental data either by 
our measurement procedure or by direct subtraction. In 
all further discussions we concentrate on NMOR signals 
only. 

Let us first study the modification of the polarization 
rotation by measuring the rotation rate 773 (-B = 0) for 
different degree of light ellipticity. We find the rotation 
rate by modulating the magnetic field by a small amount 
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FIG. 7: The normalized slope of the nonlinear magneto-optic 
rotation as a function of the ellipticity of the incident light. 
Experimental data are shown for opposite values of elliptic- 
ity and two different values of laser power: P=2 mW (solid 
up triangles for positive ellipticity and solid down triangles 
for negative ellipticity) and P=l mW (hollow up triangles for 
positive ellipticity and hollow down triangles for negative el- 
lipticity). The results of the numerical simulations for the 
case of 2 mW laser power are shown by a solid line. Abso- 
lute values of the nonlinear Faraday rotation for the linear 
polarization were d<j>/dB(B = 0) = 4.5 rad/G and 6 rad/G 
for P=2 mW and P=l mW respectively. 



FIG. 8: The normalized slope of nonlinear magneto-optic ro- 
tation as a function of the ellipticity of the incident light for 
the A scheme (transition F = 1 — * F' = 1) and M + A scheme 
(transition F = 1 — » F' = 1). Dotted lines are to guide the 
eyes. Input laser power is P=2 mW, the atomic densities are 
chosen to provide 85% absorption on each transition. The 
absolute value of the nonlinear Faraday rotation of linear po- 
larization were d(f>/dB(B = 0) = 1.8 rad/G and 4.5 rad/G for 
the F = 1, 2 — > F' = 1 transitions respectively. 



and dividing the difference of two rotation signals cor- 
responding to the small variation of the magnetic field 
by the magnitude of this variation. This way we detect 
only the rotation which depends on the external magnetic 
field. 

The rotation rate as a function light ellipticity is shown 
in Fig. 0. We observe a polarization rotation enhance- 
ment as predicted theoretically. At the same time, the 
experimental data cannot be fitted using Eq. ( ff8| ) be- 
cause of the Doppler broadening of the transition and 
the ac-Stark of the magnetic sublevels. However, an ex- 
act numerical simulation based on steady state solution 
of Maxwell-Bloch equations for the F — 2 — > F' = 1 
transition, which takes into account these effects, is in 
excellent agreement with the experimental data. 

It is also possible to verify that there is no polarization 
rotation enhancement in an isolated A system. To do that 
we tune the laser to the F = 1 — > F' = 1 transition of the 
87 Rb D\ line. In this case, the ground-state coherence is 
formed by only one A link. The relative rotation rate for 
F = 1, 2 — * F' = 1 transitions are presented in Fig. [|. 
Although there is a slight dependence of the rotation an- 
gle on the light ellipticity for F = 1 — > F 1 = 1 transition, 
this deterioration may be determined by Doppler broad- 
ening, ac-Stark shifts, etc. 

It is important to point out that even though the theo- 
retical expression for the relative rotation rate (Eq. (f7S|)) 
does not fit the experimental data precisely, it correctly 



predict some of the rotation properties. For example, our 
experiments confirm that the relative rotation rate does 
not depend on the sign of the ellipticity (Fig. ^). If we 
vary the total laser power or the coherence decay rate 
7o (by varying the laser beam diameter), the absolute 
value of the rotation changes according to Eq. (|77|); its 
dependence on the light ellipticity is the same within the 
experimental uncertainty (Figs. and |^). 

All previous data were obtained for optically thin Rb 
vapor (transmission I ut/^in — 0.85). The dependence 
of the relative rotation rate on the ellipticity for higher 
atomic densities is shown in Fig. |l(| It is easy to see that 
for nearly circular polarization the rotation decreases as 
atomic density is increased. This may be caused by op- 
tical pumping to the other ground state hyperfine levels, 
as well as by the destruction of atomic coherence by ra- 
diation trapping |i6|, |67|. 

The precise value of the output ellipticity of the laser 
polarization is required for a ccur ate polarization rota- 
tion measurements (see Eq. ( |l00|) ). The experimental 
observations demonstrate that for optically thin media 
the ellipticity of the light does not noticeably change due 
to propagation effects if the magnetic field is small. As 
the atomic density increases, however, the ellipticity in- 
creases (Fig. [To|b) . Although this change is relatively 
small (< 15%), the associated error in the calculated ro- 
tation is very significant. 
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FIG. 9: The normalized slope of nonlinear magneto-optic ro- 
tation as a function of the ellipticity of the incident light 
for two different beam diameters: d — 2 mm (circles) and 
d = 10 mm (diamonds). In both cases the laser power is kept 
at 2 mW. Absolute values of the nonlinear Faraday rotation 
for the linear polarization were d(j>/dB(B = 0) = 4.5 rad/G 
and 30 rad/G respectively. 



C. Polarization rotation of elliptically polarized 
light for large magnetic fields 

Now let us consider the case of large magnetic fields. If 
the laser frequency is swept across the atomic transition, 
the following effects contribute to the polarization rota- 
tion: nonlinear Faraday rotation due to the A-scheme 
(experimentally measured for linear polarization), self- 
rotation of elliptical polarization due to ac-Stark shifts, 
and the magneto-optic rotation of elliptical polarization 
due to M-scheme induced coherence. All these compo- 
nents are shown on Fig. [ll]. It is important to point out 
that this "new" rotation is comparable with the polar- 
ization rotation for the linear polarization and the self- 
rotation, even though this effect is due to higher order 
nonlinearity. This proves the effectiveness of the M level 
scheme for the enhancement of nonlinear susceptibility 
in atomic media. 

The magnetic field dependence of the rotation due to 
the "M-scheme" ground-state coherence reveals a very 
peculiar behavior. When the rotation is independent of 
the sign of the ellipticity in the vicinity of zero magnetic 
field (as it was demonstrated earlier), then for larger mag- 
netic fields the rotation becomes asymmetric with respect 
to both magnetic field and ellipticity. To invert the sign 
of the rotation, both the ellipticity and the magnetic field 
should change their signs (Fig. |12|a). The ellipticity of 
the outgoing light also changes with the magnetic field; 
although it is equal to the initial ellipticity for small mag- 
netic fields (at least for optically thin samples), it grows 
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FIG. 10: (a) The normalized slope of nonlinear magneto-optic 
rotation as a function of the ellipticity of the incident light for 
various atomic densities. Laser power is 2 mW, beam diame- 
ter d = 2 mm. Inset: Absolute value of the nonlinear Faraday 
rotation of linear polarization as a function of atomic density, 
(b) The output ellipticity e as a function of the ellipticity of 
the incident light for various atomic densities. Dotted line 
is for unchanged ellipticity. Inset: Transmission I Q ut/Iin of 
linear polarization as a function of atomic density. 



symmetrically when the magnetic field becomes larger 
(Fig. |l2|a). These changes must to be taken into account 
when the polarization rotation angle is measured. 



D. NMOR for atoms with higher angular 
momentum 

As discussed in Sec. IV, higher orders of nonlinear 
susceptibility may be enhanced in multi-A systems. In 
practice this means that atoms with larger ground state 
angular momentum are required. The most convenient 
candidate for the study of the higher orders of Zccman 
coherence is the 85 Rb isotope, since the same laser may 
be used as for our previous study of 87 Rb. In our exper- 
iments we use the 5S 5/2 F = 3 -► 5F 3/2 F' = 2 of 85 Rb. 
The interaction scheme of elliptically polarized light with 
this transition consists of an M scheme and a triple-A 
scheme. 

The relative rotation rate for this transition as a func- 
tion of the light ellipticity is shown in Fig. [fj|. The po- 
larization rotation enhancement, observed in this case 
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FIG. 11: The polarization rotation angle as a function of 
laser detuning for ellipticity e = 25° and magnetic field 
B — 0.35 G. The components of the rotation due to vari- 
ous processes are also shown. Zero detuning corresponds to 
the F = 2 — + F' = 1 transition. The small peak on the right 
is due to contamination of the cell with 85 Rb isotope. 
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FIG. 13: The normalized slope of nonlinear magneto-optic 
rotation as a function of the ellipticity of the incident light 
for the F = 3 — > F' — 2 transition of 85 Rb (diamonds), and 
for F = 2 — > F' = 1 transition of 8r Rb (circles). Input laser 
power is P=2 mW, the atomic densities are chosen to pro- 
vide 85% absorption on each transition. Absolute values of 
the nonlinear Faraday rotation for linear polarization were 
d(j}/dB(B = 0) = 2.9 rad/G and 4.5 rad/G respectively. In- 
set: the theoretical dependences for naturally broadened Rb 
isotopes, from Eqs. (Q) and (jsi|). 
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FIG. 12: (a) The polarization rotation angle as a function 
of magnetic field for opposite values of ellipticity. (b) The 
ellipticity of the transmitted light as a function of magnetic 
field. Initial ellipticity is shown as a dashed line. 



cited states, which may significantly change the proper- 
ties of the system. That is why it would be very interest- 
ing to measure the rotation due to high order coherence, 
discussed above, in a cloud of cold atoms. In this case 
we expect to see a much stronger effect (Eq. @)), since 
all problems caused by the overlapping transitions due to 
the motion of the atoms would be eliminated in cold gas. 

The spectral dependence of the rotation of the ellipti- 
cal polarization on laser frequency for the case of large 
magnetic field is shown in Fig. |l4|. Similarly to the 87 Rb, 
the high-order Zceman coherence significantly modifies 
the rotation spectra, and the contribution of the nonlin- 
ear rotation is comparable with the rotation of the linear 
polarization and self-rotation. 

One can see additional sub-Doppler structure on top of 
the rotation resonances. These peaks appears due to the 
retro-reflection of the laser beam inside the atomic cell. 
This additional beam interacts with atoms and causes 
the redistribution of the atomic population similar to 
Dopplcr-free saturation spectroscopy. 



is noticeably smaller than for 87 Rb. The reason for 
this may be the smaller hyperfine splitting of the ex- 
cited state (362 MHz vs 812 MHz for 87 Rb), which 
is completely overlapped by the Doppler broadening 
(^Doppier ~ 500 MHz). This overlap results in efficient 
"mixing" of the coherences induced through different ex- 



VIII. CONCLUSION 

We have studied the nonlinear magneto-optic rota- 
tion of elliptically polarized light interacting with various 
transitions of rubidium atoms. We have shown that this 
rotation can be described by means of A, M , and higher 
chain A schemes. For the simple three-level A scheme, the 
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FIG. 14: The polarization rotation angle in 85 Rb as a func- 
tion of laser detuning for ellipticity e = 25° and magnetic field 
B = 0.35G. The components of the rotation due to various 
processes are also shown. Zero detuning corresponds to the 
cross-resonance F — 3 — > F' — 2.3 transition. The distortions 
of the resonances are due to reflected light beam. 



moment. 

Since the modification of NMOR is associated with 
an enhancement of nonlinear atomic susceptibility, we 
have analyzed the effectiveness of this process by com- 
paring the nonlinear susceptibility for M and N inter- 
action schemes. We have demonstrated that although 
the enhancements of nonlinearity in these schemes are 
caused by different mechanisms, they exhibit the same 
absorptive and refractive nonlinearity magnitudes. We 
have also shown that the generalized M scheme may be 
used to create resonantly enhanced nonlinear suscepti- 
bility of any given order, similarly to the generalized N 
scheme j4j| . We have discussed the possible implementa- 
tion of the generalized M scheme for quantum computer 
algorithms. 

To verify our theoretical calculations, we have stud- 
ied the polarization rotation of elliptically polarized laser 
light propagating through Rb vapor. The M interac- 
tion scheme is realized on the F = 2 — ► F' = 1 transi- 
tion of 87 Rb, and the triple- A scheme is observed on the 
F = 3 — > F' = 2 transition of 87 Rb. Although the exper- 
imental points cannot be fit perfectly by the theoretical 
expressions (Eqs. j32] ) and d7j])), the basic properties of 
the new rotation are confirmed. 



rotation does not depend on the light ellipticity. For more 
complicated systems, the multi-photon processes are re- 
sponsible for the creation of high-order ground-state co- 
herence resulting in a new type of ellipticity-dependent 
nonlinear magneto-optical rotation. We have derived 
simple analytical expressions for this rotation for the M 
interaction scheme (Eq. (ff^)) and we showed that this 
effect results from the coherently induced hexadecapole 
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